Evidence for granularity, anisotropy and lattice distortions in cuprate 
superconductors and their implications 
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Granularity, anisotropy, local lattice distortions and their dependence on dopant concentration 
appear to be present in all cuprate superconductors, interwoven with the microscopic mechanisms 
responsible for superconductivity. Here we review anisotropy and penetration depth measurements 
to reassess the evidence for granularity, as revealed by the notorious rounded phase transition, the 
evidence for the three dimensional nature of superconductivity, uncovered by the doping depen- 
dence of transition temperature and anisotropy, and to reassess the relevance of the electron-lattice 
coupling, emerging from the oxygen isotope effects. 
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Establishing and understanding the phase diagram of cuprate superconductors in the temperature - dopant con- 
centration plane is one of the major challenges in condensed matter physics. Superconductivity is derived from the 
insulating and antifcrromagnctic parent compounds by partial substitution of ions or by adding or removing oxygen. 
For instance La2Cu04 can be doped either by alkaline earth ions or oxygen to exhibit superconductivity. The empir- 
ical phase diagram of La 2 _ a; Sr a .Cu04 0, HEM IS IE USE! depicted in FigQJi shows that after passing the so called 
undcrdoped limit (x u ~ 0.05), T c reaches its maximum value T c (x m ) at x m k, 0.16. With further increase of x, T c 
decreases and finally vanishes in the overdoped limit x D rs 0.27. This phase transition line is thought to be a generic 
property of cuprate superconductors [lfl | and is well described by the empirical relation 



T c (x)=T c (x m ) \ 1-2 ( 1 



(i) 



proposed by Presland et al.\v^. Approaching the endpoints along the x-axis, I^-zSr^CuC^ undergoes at zero 
temperature doping tuned quantum phase transitions. As their nature is concerned, resistivity measurements reveal 
a quantum superconductor to insulator (QSI) transition in the undcrdoped limit |l2l ll3lll4l lid lid Il7lll8| and in the 
overdoped limit a quantum superconductor to normal state (QSN) transition fl2l]. 

Another essential experimental fact is the doping dependence of the anisotropy. In tetragonal cuprates it is defined 
as the ratio 7 = £ a b/£c of the correlation lengths parallel (£ a &) and perpendicular (£ c ) to Cu02 layers (a6-planes). In 
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FIG. 1: (a) Variation of T c for I^-zSr^CuCU. Experimental data taken from Q, Q H H S H H @ . The solid line is Eq.JTJ 
with T c (x m ) = 39K. (b) jt versus x for I^-zSr^CuO^ The squares are the experimental data for 7t c QQQIEQ an d the 
triangles for 7t=o SB- The solid curve and dashed lines are Eq.@ with Jt c ,o = 2 and 7t=o,o = 1.63. 
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FIG. 2: Tc(x)jT c (x m ) versus jt (x m ) / Jt (x) for various cuprate families: I^-xSr^CuCU (•, T c (x m ) = 37K, 7t c (x m ) = 
20) 0,011110 , (O, T c (x m ) = 37K, 7 T=oOm) = 14.9) [|L HgBa 2 Cu0 4+ 4 (A, T c {x m ) = 95.6K, 7 t c (a5 m ) = 27) 
[13, Bi 2 Sr 2 CaCu20 8+ «5 (*, T c (x m ) = 84.2K, 7T[ (i ra ) = 133) El, YBa 2 Cu 3 7 -a (♦, T c (x m ) = 92.9K, 7Tc (x ro ) = 8) 
123, YBa 2 (Cui- H Fe H ) 3 07-i (■, T c (x m ) = 92.5K, 7Tc (x m ) = 9)0, Yi_ H Pr y Ba 2 Cu 3 7 -i (▼, T c (x m ) = 91K, 7Tc (sr TO ) = 
9.3)0, BiSraCai-^PryCuaOg (<*, T c (x m ) = 85.4K, 7t =o (x m ) = 94.3)0 and YBa 2 (Cui_ H Zn H ) 3 7 - 4 (►, T c (x m ) = 92.5K, 
7t=o (x m ) = 91 I28H . The solid and dashed curves are Eq.©, marking the flow from the maximum T c to QSI and QSN criticality, 
respectively. 



the superconducting state it can also be expressed as the ratio 7 = A c /A a b of the London penetration depths due to 
supercurrents flowing perpendicular (A c ) and parallel (\ a b ) to the afr-plancs. Approaching a nonsuperconductor to 
superconductor transition £ diverges, while in a superconductor to nonsuperconductor transition A tends to infinity. 
In both cases, however, 7 remains finite as long as the system exhibits anisotropic but genuine 3D behavior. There 
are two limiting cases: 7=1 characterizes isotropic 3D- and 7 = 00 2D-critical behavior. An instructive model 
where 7 can be varied continuously is the anisotropic 2D Ising model[l9j. When the coupling in the y direction goes 
to zero, 7 = £, x /£,y becomes infinite, the model reduces to the ID case, and T c vanishes. In the Ginzburg-Landau 
description of layered superconductors the anisotropy is related to the interlayer coupling. The weaker this coupling 
is, the larger 7 is. The limit 7 = 00 is attained when the bulk superconductor corresponds to a stack of independent 
slabs of thickness d s . With respect to experimental work, a considerable amount of data is available on the chemical 
composition dependence of 7. At T c it can be inferred from resistivity (7 = £, a b/£,c = yPab/Pc) and magnetic torque 
measurements, while in the superconducting state it follows from magnetic torque and penetration depth (7 = A c /A a b) 
data. In Fig. [flp we displayed the doping dependence of I/7T evaluated at T c (jt c ) and T = ( 7 r=o)- As the dopant 
concentration is reduced, 727 and 7t=o increase systematically, and tend to diverge in the underdoped limit. Thus the 
temperature range where superconductivity occurs shrinks in the underdoped regime with increasing anisotropy. This 
competition between anisotropy and superconductivity raises serious doubts whether 2D mechanisms and models, 
corresponding to the limit 7-7; = 00, can explain the essential observations of superconductivity in the cuprates. From 
Fig. 2]d it is also seen that 7t (x) is well described by 

1t{x) = — , (2) 

X x u 

where 7t,o is the quantum critical amplitude. Having also other cuprate families in mind, it is convenient to express 
the dopant concentration in terms of T c . From Eqs.jp) andQ we obtain the correlation between T c and jt- 

T c {x) ( 7t (x m ) \ 2 7T,0 ,„s 
= 1 - — r^T ~ 1 ' 7t \x m ) = (3) 



T c (x m ) \ 7t (x) 

Provided that this empirical correlation is not merely an artefact of L^-zSr^CuO^, it gives a universal perspective 
on the interplay of anisotropy and superconductivity, among the families of cuprates, characterized by T c (x m ) and 
jT(x m )- For this reason it is essential to explore its generic validity. In practice, however, there are only a few 
additional compounds, including Hgl^CuO^+allil, for which the dopant concentration can be varied continuously 
throughout the entire doping range. It is well established, however, that the substitution of magnetic and nonmagnetic 
impurities depresses T c of cuprate superconductors very effectivelv|2ll l22j| . To compare the doping and substitution 
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FIG. 3: T c versus (T = 0) for heavily underdoped YBa 2 Cu 3 7 -6 single crystals, taken from Hosseini|35ll . The solid line 



is T c = 170 (l/Ac (T = 0)) 1/3 and indicates the consistency with the 2D-QSI scaling relation T c oc (1/Aj? (0)) 



z/(z+2) 



and z = 1. 



driven variations of the anisotropy, we depicted in Fig. [5]the plot T c (x) /T c (x m ) versus jt (x m ) / 7t (x) for a variety 
of cuprate families. The collapse of the data on the parabola, which is the empirical relation J3J), reveals that this 
scaling form appears to be universal. Thus, given a family of cuprate superconductors, characterized by T c (x m ) and 
7t {xm)-, it gives a universal perspective on the interplay between anisotropy and superconductivity. 
Close to 2D-QSI criticality various properties are not independent but related bv |13L UA . Ha. Ha ITa Il8l| 



T r = 



HR-2 



l6^k B Xl b (0) 



7 T 



(4) 



where k B is the Boltzmann constant, and $o the elementary flux quantum. A Q {, (0) is the zero temperature in-plane 
penetration depth, z is the dynamic critical exponent, d s the thickness of the sheets, and v the correlation length 
critical exponent of the 2D-QSI transitions. S measures the distance from the critical point along the x axis (see Fig^i, 
and i?2 is a universal number. Since T c oc d s /\ 2 ab (0) cx rig , where is the aerial superfluid density, is a characteristic 
2D property, it also applies to the onset of superfluidity in 4 He films adsorbed on disordered substrates, where it is 
well confirmed]^]. A great deal of experimental work has also been done in cupratcs on the so called Uemura plot, 
revealing an empirical correlation between T c and d s /X^ b (0)[3(|- Approaching 2D-QSI criticality, the data of a given 
family tends to fall on a straight line, consistent with Eq.Q. Differences in the slope reflect the family dependent 
value of <i s , the thickness of the sheets, becoming independent in the 2D limit [l4L lift 1171 1 181. The relevance of d s 
was also confirmed in terms of the relationship between the isotope effect on T c and 1/ \ 2 ab [ijj |3J . Moreover, together 
with the scaling form (0} the empirical relation JTJ implies 2D-QSI and 3D-QSN transitions with z = 1, while the 
empirical relation for the anisotropy (Eqs.(0) and Q), require v = 1 at the 2D-QSI critical point. Thus, the empirical 
correlations point to a 2D-QSI transition with z = 1 and V = 1. These estimates coincide with the theoretical prediction 
for a 2D disordered bosonic system with long-range Coulomb interactions, where 2=1 and V ~ 1 j32l I33I l34l| . Here 
the loss of superfluidity is due to the localization of the pairs, which ultimately drives the transition. From the scaling 
relation it is seen that measurements of the out of plane penetration depth of sufficiently underdoped systems 

allow to estimate the dynamic critical exponent z directly, in terms of T f oc (l/A^ (0)) z ^ z+2 \ which follows from 
Eq.(@J with j T = X c (0) /A ah (0). In Fig^Jwe displayed the data of Hosseini[35| for heavily underdoped YBa 2 Cu 3 07_,5 



The solid line is T c = 170 
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1/3 

1/A^ (T = 0)) and uncovers the consistency with the 2D-QSI scaling 



single crystals. 

relation T c oc (l/A^ (0))*' with z = 1. 

We have seen that the doping tuned flow to the 2D-QSI critical point is associated with a depression of T c and 
an enhancement of 7t- It implies that whenever a QSI transition is approached, a no vanishing T c is inevitably 
associated with an anisotropic but 3D condensation mechanism, because 7t is finite for T c > (see Figs^ 3 and[2j. 
This represents a serious problem for 2D models[3r| as candidates to explain superconductivity in the cuprates, and 
serves as a constraint on future work toward a complete understanding. Note that the vast majority of theoretical 
models focus on a single Cu-0 plane, i.e., on the limit of zero intracell and intercell c-axis coupling. 

Since Eq.|@J is universal, it also implies that the changes AT C , Ad s and A (l/A^ b {T = 0)), induced by pressure or 
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FIG. 4: Data for the oxygen isotope effect in underdoped La 2 - 3: Sr 3: Cu 04( 0: x=0.15|2l, A:x=0.08, 0.086 
Yi-.Pr^BaaCugOy-A- (•: x=0, 0.2, 0.3, 0.4)111 |H EJ and YBa2Cu 3 07-,5 (T 0, M^) in terms of A (\ ab (0)) / X ab {0) 
versus -AT C /T C . The solid line indicates the flow to 2D-QSI criticality and provides with Eq.@ an estimate for the oxygen 
isotope effect on d s , namely Ad 3 /d s = 3.3(4)%. 



isotope exchange are not independent, but related by 

AT C _ Ad s A (1/A*„ (0)) _ Adg A (A afc (0)) 

T c d s + (1/A* 6 (0)) d s X ab (0) • 1 ' 

In particular, for the oxygen isotope effect ( 16 vs. 18 0) of a physical quantity X the relative isotope shift is defined 
as AX/X = ( 18 X - 16 XV™ X. In FigEJwe show the data for the oxygen isotope effect in La9-.Sr-.CuO/i [37l 13^ . 
Yi_ 2: Pr 2: Ba2Cu307_5j3l,l3^,E3 an d YBa2Cu3O7_5|3l,0|, extending from the underdoped to the optimally doped 
regime, in terms of A (A a b (0)) / A a b (0) versus AT C /T C . It is evident that there is a correlation between the isotope effect 
on T c and X ab (0) which appears to be universal for all cuprate families. Indeed, the solid line indicates the flow to the 
2D-QSI transition and provides with Eq.JSJ) an estimate for the oxygen isotope effect on d s , namely Ad s /d s = 3.3(4)%. 
Approaching optimum doping, this contribution renders the isotope effect on T c considerably smaller than that on 
X a b (0). As shown in Fig|SJ even in nearly optimally doped YP^CuaOy-s, where AT C /T C = —0.26(5)%, a substantial 
isotope effect on the in-plane penetration depth, AX ab (0) / X ab (0) = —2.8(1.0)%, has been established by direct 
observation, using the novel low-energy muon-spin rotation technique j^jj. Note that these findings have been obtained 
using various experimental techniques on powders, thin films and single crystals. 

Since upon oxygen isotope exchange the lattice parameters remain essentially unaffected 0, El- the substantial 
isotope effect on the in-plane penetration depth uncovers the coupling between local lattice distortions and superfluid- 
ity and the failure of the Migdal-Eliashberg (ME) treatment of the electron-phonon interaction, predicting, 1/A 2 (0), 
to be independent of the ionic masses 01- Evidence for this coupling emerges from the oxygen isotope effect on d s , 
the thickness of the superconducting sheets, upon isotope exchange, while the lattice parameters remain unaffected. 
Indeed, the relative shift, Ad s /d s ~ 3.3(4)%, apparent in Fig|31 implies local distortions of oxygen degrees of freedom, 
which do not modify the lattice parameters, and are coupled to the superfluid. 

Further evidence for this coupling emerges from the isotope effect on the granularity of the cuprates |45l l46j | . Recently, 
it has been shown that the notorious rounding of the superconductor to normal state transition is fully consistent 
with a finite size effect, revealing that bulk cuprate superconductors break into nearly homogeneous superconducting 
grains of rather unique extent [lj, 1481. E ven evidence for their surface and edge contributions to specific 

heat and penetration depth has been established[49j. A characteristic feature of a finite size effect in the temperature 
dependence of the in-plane penetration depth X ab is the occurrence of an inflection point giving rise to an extremum 
in d (X3. h (T = 0) /A?b (T)) jdT at T p . Here X ab (T p ), T p and the length L c of the grains along the c-axis are related 
byllllil III 

1 Wn 3 k B T p 
X 2 ab {T P ) *IL C ■ (b) 

Recently we explored the effect of oxygen isotope exchange in Yi-zPr^P^CuaOr-s on L c by means of in-plane 
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FIG. 5: Magnetic field penetration profiles B(z) for a 16 O substituted (closed symbols) and a 18 O substituted (open symbols) 
YBa2Cu3C>7_,5-film measured in the Meissner state at 4 K and an external field of 9.2 mT, applied parallel to the surface of 
the film. The data are shown for implantation energies 3, 6, 10, 16, 22, and 29 keV starting from the surface of the sample. 
Solid curves are best fits to B(z) = Bo cosh[(t — z)/\ a b\ co$fa{t / \ a b) ■ This is the form of the classical exponential field decay in 
the Meissner state B(z) = Bo exp(— z/A a (,), modified for a film with thickness 2t with flux penetrating from both sides. Taken 
from Khasanov et al. Ell 



penetration depth measurements ^(|. Note that the shifts arc not independent but according to Eq.© related by 

AL C = AT Pc | AXl b (T Pc ) 
Lc T Pc X 2 ab (T Pc ) 



(I) 



From the resulting estimates, listed in Table I, several observations emerge. First, L c increases systematically with 
reduced T Pc . Second, L c grows with increasing x and upon isotope exchange ( 16 0, 18 0). Third, the relative shift of 
T Pa is very small. This reflects the fact that the change of L c is essentially due to the superfluid, probed in terms of 
X 2 ab . Accordingly, AL C /L C w AX 2 ab /X 2 ab for x = 0, 0.2 and 0.3. 



X 





0.2 


0.3 


AL C /L C 


0.12(5) 


0.13(6) 


0.16(5) 


^T Pc /T Pc 


-0.000(2) 


-0.015(3) 


-0.021(5) 


AXi b (T p J/Xi b (T Pc ) 


0.11(5) 


0.15(6) 


0.15(5) 


iB T Pc (K) 


89.0(1) 


67.0(1) 


52.1(1) 


ib i c (A) 


9.7(4) 


14.2(7) 


19.5(8) 


i8 r Pr (K) 


89.0(1) 


66.0(2) 


51.0(2) 


L *L C (A) 


10.9(4) 


16.0(7) 


22.6(9) 



Tabic I: Finite size estimates for the relative changes of L c , T Pc and X 2 b (T Pc ) upon oxygen isotope exchange for 
Y 1 _,Pr ;c Ba 2 Cu30 7 _ 5 |4i. 

To appreciate the implications of these estimates, we note again that for fixed Pr concentration the lattice param- 
eters remain essentially unaffected Oil ■ Accordingly, an electronic mechanism, without coupling to local lattice 
distortions implies AL C = 0. On the contrary, a significant change of L c upon oxygen exchange requires local lattice 
distortions involving the oxygen lattice degrees of freedom and implies with Eq.Q a coupling between these distor- 
tions and the superfluid. A glance to Table I shows that the relative change of the grains along the c-axis upon oxygen 
isotope exchange is significant, ranging from 12 to 16%, while the relative change of the inflection point at T Pc , or the 
transition temperature, is an order of magnitude smaller. For this reason the significant relative change of L c at fixed 
Pr concentration is accompanied by essentially the same relative change of X 2 b , which probes the superfluid. This 
uncovers unambiguously the existence and relevance of the coupling between the superfluid and lattice distortions, 
involving the oxygen lattice degrees of freedom. Furthermore the substantial isotope effect on the in-plane penetration 
depth at T = T Pa extends the evidence for the failure of the Migdal-Eliashberg (ME) theory of the electron-phonon 
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interaction, predicting 1/A 2 to be independent of the ionic masses [44|. to finite temperature. Although the majority 
opinion on the mechanism of superconductivity in the cuprates is that it occurs via a purely electronic mechanism 
involving spin excitations, and lattice degrees of freedom are supposed to be irrelevant, the relative isotope shifts 
AL C /L C ss AA 2 h / 16 A 2 h w 14% and Ad s /d s « 3% uncover clearly the existence and relevance of the coupling between 
the superfluid and local lattice distortions. Recent site-selective oxygen isotope ( 16 0/ 18 0) effect measurements of 
the in-plane penetration depth in Y .6Pro.4Ba2Cu307_5, using the muon-spin rotation (//SR) technique show that 
the distortions arise from the oxygen sites within the CuC>2 planes (100 % within error bar) j4jj. Potential candi- 
dates are then the Cu-0 bond-stretching-type modes showing a temperature dependence, which parallels that of the 
superconductive order parameterjfifij. 

To summarize we observed remarkable consistency between the scaling properties of the experimental data for a 
variety of cuprates and those characterizing 2D-QSI transitions. The important implication there is that in cuprates 
a no vanishing transition temperature and superfluid density in the ground state are unalterably linked to a finite 
anisotropy. Furthermore, the oxygen isotope effect on the in-plane penetration depth and the spatial extent of the 
superconducting grains revealed the coupling between local lattice distortions and superfluidity, while the lattice pa- 
rameters remain essentially unaffected. These findings raise serious doubts that 2D models i neglecting granularity 
and local lattice distortions are potential candidates to explain superconductivity in cuprates. 

The author is grateful to D. Di Castro, R. Khasanov, K.A. Miiller, and J. Roos for very useful comments and 
suggestions on the subject matter. This work was partially supported by the Swiss National Science Foundation and 
the NCCR program Materials with Novel Electronic Properties (MaNEP) sponsered by the Swiss National Science 
Foundation. 
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